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Abstract 

The dissipative dynamics of a heavy quark passing through charged thermal plasmas of 
strongly coupled A/" = 4 super Yang- Mills theory is studied using AdS/CFT. We compute 
the linear response of the dilaton field to a test string in the rotating near-extremal D3 brane 
background, finding the momentum space profile of (tiF^) numerically. Our results naively 
support the wake picture discussed in hep-th/0605292, provided the rotation parameter is 
not too large. 
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1 Introduction 



According to the AdS/CFT correspondence [H 13 El El, the physics of a heavy quark passing 
through finite temperature A/" = 4 super Yang-Mills plasmas may have the dual description 
in terms of a fundamental string in the background formed by a stack of near-extremal D3- 
branes. Such a subject has been studied in a number of recent papers [3 El 13 El El UHl 
im 1121 CHI UH Uni Uni, motivated by its possible connection with jet-quenching observed at 
RHIC fTlElElEni in relativistic heavy ion collisions. This remarkable phenomenon can be 
understood as the strong energy loss of a high energy parton moving through the quark-gluon 
plasma; for some non-AdS/CFT based theoretical studies, see e.^. [HI E21 El ESI El EHl EZ| • 

In jH] in, the AdS/CFT duality was applied to the computation of the drag force on a mov- 
ing quark received by a hot A/" = 4 SYM plasma. The background metric in this dual descrip- 
tion comes from near-extremal static D3-branes, which is simply ArfS's-SchwarzschildxS'^ in 
the near-horizon limit and corresponds to a neutral quark-gluon plasma. Such a computation 
was extended to i?-charged A/" = 4 SYM plasmas [H] , where the background is taken to 
be the near-horizon geometry of rotating D3 branes, giving rise to the drag force exerted by 
the thermal plasma with a non-vanishing chemical potential. However, both of these works 
relied on the test string approximation; back-reaction of such a string on the background 
was completely neglected. More recently, the authors of J2] have computed linear responses 
of the dilaton field to a test string in the ArfS's-SchwarzschildxS'^ background, which takes 
some back-reaction effects into account. Their result is quite interesting, exhibiting a recoil 
energy scale in some ranges and a directional peak of gluons radiated from the heavy quark, 
in consistent with the wake scenarios [23 EZI of jet-quenching. With AdS/CFT, this kind 
of computations should provide us useful information on the energy loss of heavy quarks 
passing through the plasma. 

In this note we wish to extend the result of ^2] to -R-charged thermal plasmas, by studying 
the linear response of the dilaton field to a test string in the rotating near-extremal D3 brane 
background. We shall treat the neutral plasma as a special case of the charged ones, to see 
how the profile of (trF^) changes and, in particular, whether a wake form can be found in 
the presence of charges. For simplicity, we will restrict ourselves to the case where only one 
angular momentum parameter is no n- vanishing j2ni|Sni- The near horizon geometry thus 
reads 



ds'^ = f-\[-hdt^ + dx^) + f^h-^dr^ 
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with 



J HA' 
A = 1 + — cos^ 6*, 
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and 



gfy^A^a'^. For later convenience, let us introduce two symbols hi, h2 
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where hi is the event horizon. We will consider a classical string placed in the background 
(0), with one end attached to the boundary of the AdS space moving in a constant velocity 
u, and the other end attached to the horizon of the black hole at the center of the AdS space. 
Our purpose here is to compute the linear response of the dilaton in the background metric 
(H}, paying particular attention to the energy scale at which dissipation occurs as well as the 
existence of directional peaks. 

This paper is organized as follows: In section 2 we recall some relevant formulae of pilllj. 
In section 3, we derive a set of differential equations governing linear perturbations of the 
dilaton field, following ^21 closely. Section 4 then applies the relaxation method [2H] to the 
boundary value problem, and presents our numerical results with some discussions. As we 
will see, in the case of / = 0, our results agree with those given in ^2]; where a different 
numerical method was used. 



2 Test String Solutions 



Let us begin with the string configuration 



r, X' = a, X^ = vT + e(r), 6 = 9{r) 



(4) 
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The function ^(r) in (@]) is determined by 



hh 



(5) 



and 6{r) obeys a complicated equation [HIE]; which in generical (for / 7^ 0) does not allow 
constant solutions. However, there are two special cases where 6{r) can have r- independent 
solutions, provided ^ ^ or | as r ^ 00. In the first case the string is parallel to the rotation 
axis, while in the second case the string lies in the plane perpendicular to the rotation axis. In 
the present paper we will consider the first case only. This string configuration (with 6 = 0), 
though special, is nevertheless non-trivial, since it could be used to study the physical effects 
of the SYM plasma charge on the test string. 
Inserting 9 = into ^ gives 



^ L2 h 



r4 + r2/2 
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One may write down an explicit expression for ^(r) by integrating out the above equation. 
With the definition Pa |) -(|3E |) . it is easy to derive: 



e(r) 
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1 , + hi 1 _ r 
— in —^=^= — — arctan — 
hi Vr - hi h2 h2 
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Now, to compare this with jH] f2] , we set 
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so that 
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irhi r + hi 2hi r 
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h2 r - hi /i2 h2 



(9) 



Note that Eq.(jn} under the limit / —* will become identical to the string configuration 



e(r) 
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studied in |H]U21- This could be expected, since in that limiting case the metric reduces 
exactly to ArfS's-Schwarzschild. Now for generic / 7^ 0, the effective 5-dimensional background 
felt by the string parallel to the rotation axis also gets simplified. Actually, taking 6 = 



3 



and freezing all the angular degrees of freedom in (^J yields: 



/ 



h= 1 
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3 Linear Responses of the Dilaton 

We turn now to back-reaction of the test string to the background, following the same routine 
of ^2|- Consider the response of the dilaton field to the test string, where has non-trivial 
solutions determined by minimizing the action 



S 



44 



dx^V^{d(l)f 



d^ae^ 
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with 



2 r 3 



-fab = daX^d.X'^G.,, 4 = SttG, 



(12) 



To establish the equation of motion for 0, it is convenient to express the action as a single 
volume integral 



S 



1 



44 



A4_ 



d ae 



M 
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This gives rise to the linearized equation of motion describing the response of the dilaton to 
the test string: 



1 



2-Ka' 



£ae^nX^S\x^ - X^(a)) = J 



(14) 



Substituting (jlj), (jH)), into the above definition of J and calculating v^— 7 as well as \/ —G 
in terms of the metric components, one finds 



J 
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2% a' 4^ 

= = ^^/'(^ = 0) 



b{x' - X\t, r))5{x' - Xy{x' - X') 



(15) 
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Suppose, just as in [T2], that depends on and t only through the combination —vt. 
Also, notice that in our string configuration, ^ = 0. Hence we can simplify Eq.(|14j) to 



27ra 



- — b{x' -vt- i{r))5{x' - X')5{x' - X^) 



(17) 



Now, going to the momentum space 



0(t,r,x) 



(27r)3 

the above equation is transformed into 



2Txa' 



dr[f ^h{r)drMr)] - f{r) 



(19) 



here / and g are given in (fTn|l . 

It seems not possible to solve Eq. (jl9j) analytically for the full range of r, but nevertheless 
we can find explicit solutions in two asymptotical regions. The first one is in the near-horizon 
region, where r ~ /ii and therefore 
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In this limit ()19|) reduces to 



4^9,(r - hi)drMr) + Trr^kfMr) 

Lho{r — hi) 



exp 
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Now let 



F = ln^, P 
hi 



TThi 2hi hi 

— m 2 ; — arctan — , 

h2 n.2 h2 



(22) 
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we can rewrite the above differential equation as 

vkiZH 



,2 
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The general solution of Eg .(1231) thus takes the form 



z nhiL ^Y^-ivkiZH(Y+p)/4. /23\ 
4r^ ^ ^ 



kNHir) = , • . /o + C^e-'^^-"-/' + C,-e-^--^/^ (24) 

4rg 1 - ivkiZH/2 " " 

where are arbitrary constants. In this near horizon region, we require (f) depend on t and 
Y through the combination t + vkiZnY/A only, so we can set = 0. Physically this means 
that we accept the infalling solution while reject the outgoing one. 

Next we consider asymptotical solutions in region near the (AdS) boundary, where 

r ^ oo, /"2 ^ ^, /i ^ 1 (25) 



In that region, Eq.(fT^ becomes 



-OyOrMr) = — (26) 
r r 

whose general solution has the form: 

4>k,NB{r) = + + B^r-^ (27) 

The constant should be set to zero, as there are no deformations in the dilaton Lagrangian. 
Physically we will be interested in Bk, since according to AdS/CFT, this quantity is directly 
related to the vacuum expectation value of the operator (9^2 ~ trF^ coupled to the dilaton. 
Actually, using the AdS/CFT dictionary we can write the VEV as 

1 1 

{OF-2{t,x)) = — ^lim^^ooV-^fi''"''S^0 = --T-^limr^oc-j-^drCpit, x,r) (28) 

where the metric (liup has been applied. Note that the above limit does not exist due to the 
first term in Eq. (j27|) . Since such a term does not depend on k, transforming it into real space 
will lead to a delta function supported at the quark location, which should be subtracted 
j. The correspondence between {Of2) and 4>(t,x,r) in should now be understood as 
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a relation after the delta-function subtraction. We thus find: 



{O^^it^x)) = J A^e*(^^-^*)+'^-^+^'=3.35^_ (29) 

Now, as in the neutral plasma case jT2j, further subtraction is needed in order to separate 
the dissipative dynamics from the near field contributions of the quark. One expects that 
these near field contributions correspond to the string hanging straight down in AdS^. Thus, 
when f = 0, we can apply the method of [31] directly to derive 



^^'^ = Y^^^^^- (30) 



Comparing this with for w = 0, one gets 



In the case of f 7^ 0, the near field contributions are obtainable through a Lorentz boost to 
the string configuration. The result reads 



2 



nncar field " ^ 



^il-v^)kl + kl, (32) 



16 

which will be subtracted from the numerical values of Bk computed in the next section. 

4 Numerical Results and Discussions 

We are now in a position to solve the following boundary value problem numerically: 



drlf ih{r)dr(j)k{r)] - f{r) 



fc(r) = e-'"'^^'' (33a) 



Mr) l-,,fc,,,/2 + ^ ' = + ^33b) 

Mr) ^ Bkr~^ = (f)k,NBl + <i>k,NB2 (33c) 

To this end, we need to implement the boundary conditions at r ~ /ii and r ~ 00, and this 
can be done by introducing two Wronskians, W^nir) and WiyB^r), to measure the differences 
between our numerically evaluated 0fc(r) and the asymptotical solutions (j2D), (j77j) found in 
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the last section. One thus defines, following [T^ . 

WNH{r) = {(pk - <Pk,NHl)(p'k,NH2 " - <P'k,NHl)'Pk,NH2 
WNsir) = {^k - (i>k,NBl)4>'k,NB2 - i<P'k - <P'k,NBl)(i>k,NB2 (34) 

The boundary conditions can then be imposed properly by the requirements WNH{r) = at 
r ~ /ii and W]\fB{r) = at r ~ cxd. 

The numerical recipes we used in solving the boundary value problem (j33|) is the standard 
relaxation method j2HI- In programming, we changed the variable r ^ y = hi/r, so that 
the integration range becomes y G (0, 1), where y ~ 1 corresponds to the horizon and y ~ 
to the AdS^ boundary. Since the solution has oscillating behavior at y ~ 1, we divided the 
the integration region into four intervals, (1,0.9], (0.9,0.7], (0.7,0.4] and (0.4,0), and then 
divided the first, second, third, and forth intervals into 4096, 1024, 256, and 64 integration 
steps, respectively, in order to reach as high precision as possible while keeping CPU time 
in an acceptable range. We also checked that if the whole integration region is divided into 
50000 steps uniformly, one can get almost the same results and, of course, this costs more 
CPU time. 

Since Eq. ()33a|) becomes singular as r — > /ii or r — > oo, we have to set Wnh = at a point 
very close to the horizon, and set Wnb = at a large but finite value of r. The conditions 
imposed in our computations are 

WNHl=imih, = 0' WMBl=7ooh^ = 0- (35) 

We depict our numerical results in FigCJI^l where equi- value lines of Re-B^, —ImBk, 
etc., for some different values of v and /, are plotted in the momentum plane k = [Ki, K±). In 
all the plots the near field contributions ()32|) have been subtracted. We adopt the convention 
of ^21; representing values closest to zero by white regions, and representing the most positive 
values by black regions. To compare our results with those of ^21; we fixed the energy scale 
explicitly, by setting the temperature of the plasma 

to be T = l/yrGeV = 318MeV, in accordance with the choice made in ^21- Thus, the wave 
numbers k displayed here are measured by Gev/c. Apart from this, the rotation parameter 
/ is shown in units of hi, namely, when we say / = 0.5, we actually mean / = 0.5hi. 
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Figure 1: Distributions of equi-value lines of KeB{Ki, K±), —lmB[Ki,K±), \B{Ki, K±)\, 
and K±\B{Ki, K±)\ in the momentum plane {Ki,K±), for angular momentum / = 0, at 
speeds v = 0.75, 0.90, 0.95, 0.99. We have subtracted the near field contribution (jH^ from 
B{Ki,Kj_). 

Fig^shows our results of Bk = B{Ki, K±) in the special case / = 0. This corresponds to 
the neutral plasma studied in ^21- The first two columns contain plots of IleB{Ki, K±) and 
—lmB{Ki, K±) at V = 0.75, 0.90, 0.95, 0.99, while the third and forth show the corresponding 
plots of \B{Ki, K±)\ and K±\B{Ki, K^)\, respectively. Comparing these with the plots given 
in one sees that the basic features are the same, both exhibiting a directionally peaked 
structure in K±\B{Ki, K±)\ and a possible range of the recoil energy. For example, focusing 
on the third line, third column of Fig^ we find that \B{Ki, K±)\ (for v = 0.95) is peaked at 

K^^O, 3<Ki< 7.2 Gev/c, 

which indicates that the recoil energy Er is in the range 1.5 < Er < 3.6 Gev, less than the 
value 

Ef = —— T = 6.2 Gev (for v = 0.95) 

\ — 

predicted in the corresponding free field theory by a factor of a few. This agrees perfectly 
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Figure 2: Distributions of equi-value lines of \B{Ki, K±)\ and K±\B{Ki, K±)\ in the plane 
{Ki,K±), for angular momentum / = 0.5 and 1, at speeds v = 0.75,0.90,0.95,0.99. We 
have subtracted the near field contribution ()32|) from B[Ki, K±). The left eight ones are 
plots of \B{Ki, K±)\ and K±\B{Ki, K±)\ for / = 0.5hi, and the right eight ones are the 
corresponding plots for / = l.Ohi. 

with the numerical result of [T^ . 

Now we turn to the case / 7^ 0. When / = 0.5 and 1, our results for B{Ki, K±) (at speeds 
V = 0.75, 0.90, 0.95, 0.99) are displayed in Fig|21 The first and third columns give the plots 
of \B{Ki, K^) \ for / = 0.5 and / = 1, respectively. From these plots, one sees that there may 
exist a finite range of K± in which \B{Ki, K±)\ is peaked, just as in the / = case. Taking 
the plot for 1 = 1 at = 0.95 as an example (the one placed at the third line, third column 
of Fig|2))), we find that \B{Ki, K^)\ has a peak within 

K^^O, 2<Ki< 7.6 Gev/c, 

so that the recoil energy is roughly in the range 1 < Ej. < 3.8 Gev, which is slightly larger 
than the result found in the / = case. That the range of E^ (in particular its upper bound 
value) becomes larger when / increases seems to be a generic phenomenon in our numerical 
computations. A possible implication is that the more is the charge carried by the plasma. 



10 



the more energy of the quark would be dissipated away by radiation of gluons. This should 
not be taken too seriously, however, since our numerical results may contain more errors as 
/ becomes larger; see the discussion for / = 2 below. 

The forth column of Fig|21 shows plots of K^\B{Ki, K±)\ for / = 1, at speeds v = 0.75, 
0.90, 0.95, and 0.99. Again, we can clearly see a directionally peaked structure in these plots, 
much resembling the / = case. We have also displayed plots of Kj_\B{Ki, K^)\ for / = 0.5 
in the second column of Fig|21 finding that they look quite similar to (and in fact, they are 
intermediate between) those in the / = 0, 1 cases. Thus, at least for / < 1, our results suggest 
that the wake picture described in [T2] may also apply to i?-charged plasmas. 




Figure 3: Distributions of equi-value lines of \B{Ki, K^)\ and K±\B{Ki, K^)\ in the plane 
{Ki,K±), for / = 2, at speeds v = 0.75,0.90,0.95,0.99. The near field contribution §^ 
has been subtracted from B{Ki, Kj_). The left eight plots are computed using the boundary 
conditions (j35|) . while the right eight ones are the corresponding plots with the new boundary 
conditions WNH\r=imihi = 0, WNB\T=soohi = imposed. 

What would happen if I becomes larger? For a comparison, the plots of \B{Ki, Kj_) \ and 
K±\B{Ki, K±) \ for / = 2 are displayed in the first and second columns of Fig|21 The pattern 
of such plots appears now to be somehow different from what we have seen in the case of / < 1. 
In particular, at relatively low speeds, the directionally peaked structure in K±\B{Ki, K±)\ 
seems no longer obvious. We should mention that in our numerical calculations, plot patterns 
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for / being larger will depend more sensitively on how we impose the boundary conditions. 
For instance, if we modify the second condition in (|35|) by setting W^b = at r = 800/ii 
instead of TOO/ii, then the plots for / = 2 will have a new pattern, as shown in the third and 
forth columns of Figl^l which is different from the previous one. Similar changes in plots 
are also observed at / = 0, 0.5 and 1, but they are far less remarkable, leaving the pattern 
qualitatively the same. Hence, our numerical results for / = 2 are not quite reliable, which 
may contain more errors than those in the / < 1 case. 

Following we can determine the opening angle 6 between the velocity of the heavy 
quark and the directional peak found in K^B{Ki, K^). The results are summarized in the 
following table: 



e 


V = 0.75 


V = 0.90 


V = 0.95 


V = 0.99 


I 


= 


0.59 


0.42 


0.28 


0.17 


I 


= 0.5 


0.57 


0.39 


0.27 


0.16 


I 


= 1 


0.55 


0.37 


0.26 


0.15 



We thus find a strong dependence of 6 on v in both the neutral and charged cases. For fixed 
V, however, the opening angle depends rather weakly on / (at least in the region / < 1; as 
we mentioned, our numerical results for I = 2 are not quite trustable), in a monotonically 
decreasing way. 

In conclusions, we have studied the dissipative dynamics of a heavy quark passing through 
charged M = 4 SYM plasmas. Neutral plasmas were treated as a special case of the charged 
ones, where we reproduced the main results of ^21 using a different numerical method. Our 
results for / < hi naively support the wake picture, but they are inconclusive for / ~ 2hi due 
to numerical errors. 
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